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Systems with itinerant fermions close to a zero temperature quantum phase 
transition like the high temperature superconductors exhibit unusual non- 
Fermi liquid properties. The interaction of the long-range and low-energy 
fluctuations of the incipient order with the fermions modify the dynami- 
cal properties of the fermions strongly by inducing effective long-range in- 
teractions. Close to the transition the interaction of the order parameter 
fluctuations becomes important. In this paper we discuss the physics of the 
non-Gaussian order parameter fluctuations on the electronic spectrum and 
illustrate their effect by considering the charge-density-wave transition and 
a two dimensional superconductor. 
PACS numbers: 05.70 Ln, 05.70 Jk, 64. 



1. INTRODUCTION 

Fermionic systems which undergo a zero temperature phase transition 
as a function of some experimental parameter exhibit unusal physical proper- 
ties. The phase transition which are driven by quantum rather than thermal 
fluctuations have been studied intensively theoretically and experimentally 
the last few years. While significant progress in the understanding of low 
dimensional spin systems could be achieved the physics of itinerant electrons 
coupled to the low energy collective fluctuations is far less understood. 

One of the problems studied most intensively in this context is a itiner- 
ant electron magnetism. The interplay between superconductivity and spin 

udied in nearly 
used the term 



fluctuations in nearly ferromagnetic metals has been first t 
ferromagnetic metals. Berk and Schriefferi and DoniachI 
"paramagnons" to describe the long-lived and long-range spin fluctuations 
in a nearly ferromagnetic metal. The exchange of long-lived and long-range 
ferromagnetic spin fluctuations has been argued to mediate an unusual type 
of superfluidity in ^He. The nonlinear feedback of spin-fluctuations has been 



used sucessfully to explain the stability of the " Anderson-Brinkmann-Morel" 
superfluid phase. However by now it has bectinie clear that ^He closer to 
localization than to a ferromagnetic transitioned. In recent years it has been 
proposed that in the high temperature superconductors with their magnetic 
parent compounds the exchange of anti ferromagnetic spin fluctuations may 
lead to d-wave superconductivity taliJ. Some experimental features like the 
the unusual relaxation rates in the NMR experiments have beenunderstood 
using a phenomenologicalical approach to the spin fluctuationgiil. Various 
re-summation schemes, like the fluctuation-exchange approximation (FLEX) 
and other more complicated schemes like Parquet re-summationu, have been 
used to calculate properties in the "nearly antiferromagnetic phase" and to 
compare it to more microscopic approaches. Already in 1976 Hertz showed 
how to generalize the concepts of classical phase transitions to a zero tem- 
perature quantum phase transition for itinerant electron^. While thermo- 
dynamic properties of the electronic system can be calculated close to the 
phase transition the spectral properties of the electrons have been calcu- 
lated only in approaches basically relying on the assumption of Gaussian 
order parameter fluctuation^. Very little is known about the physics of 
the electronic degrees of freedom in the non-Gaussian regime. The purpose 
of the paper is to draw attention to the fact that the non-Gaussian order 
parameter fluctuations drastically modify the electronic properties. This is 
particularly relevant for the pseudo-gap regime. The paper is organized as 
fowllows. First we briefly review the theory of quantum phase transitions 
by Hertz and Millis. We then explain the physics for of the non-Gaussian 
regime using a one-dimensional example. Finally we apply how to general- 
ize the one-dimensional example to the problem of phase fluctuations in a 
two-dimensional superconductor. 



2. QUANTUM PHASE TRANSITIONS 

The partition , Z, function of an interacting electron gas with incipient 
order can be expressed as a functional integral over a field $ which represents 
the order-parameter density of the incipient order, and the electronic degrees 
of freedom, iJJ^: 

Z = Z^ j Vm I V[ij] exp(-5[V^, $]) (1) 
The last integal defines an effective action, 5'e//[<&], for the bosonic field 



(2) 



The electronic degrees of freedom can be integrated out, since they only 
appear quadratic in the action. An external field hn{q) can be coupled to 
the field $ via 

si}f=PY.K{q)-^-n{-q) (3) 

Integrating out the electronic degrees of freedom and expanding in the field 
$ the action S[^] can be written as 

Seff m = sf^f + sf]f m + sfjf m + ... (4) 

as shown in Fig. ||. The zig-zag line represents the field <I> and the solid line 
the electron propagator. 
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Fig. 1. Expansion of the partion function Z in terms of the Hubbard- 
Stratanovich fields. 



The second order term is given by: S'(^)[^>] = /3 X;i,2 -f'^^^Hl. 2)$(1)^>(2) 
where 1 represents (gi, zwi, ai) and so on. The sum over the bosonic Mat- 
subara frequencies iuj is multiplied with the inverse temperature (3. K'^'^\1, 2) 
can be calculated explicitely K^'^\l, 2) = 5^^i+ui2^qi+q25a,a'U (1 — Uxoiliii^i)) 
where xo(9ii ^"^i) is the electronic susceptibility and U is the coupling strength 
of the electrons to ^. 



€p+g - ep 



(5) 



where /(e^) is the Fermi distribution function. In a paramagnon type of 
theory only K^^^ is considered and the interaction of the the bosonic fields is 
neglected. In the self-consistent renormalized (SCR) approach by Moriya 
the interaction of the bosonic modes is taken into account by renormalizing 



K'^'^'^ . The SCR therefore an effective gaussian theory. The effects we are 
considering here are due to the interaction of the bosonic modes which is 
given by the next term in the series for Seff- 

The fourth order term obviously depends on four bosonic fields and is 
given by: S^^\^ = /3 X]i,2,3,4 2, 3, 4)$(1)$(2)$(3)$(4). If the Fermi 

surface is not nested and the ordering wave vector of the incipient order is 
not an extremal spanning vector of the Fermi surface then the complicated 
fourth order term might be approximated by a constant. The effective action 
including the fourth order term was first studied by Hertz 11 and subsequently 
by MillisH using renormalization group. The resulting phase diagram for a 
two dimensional system is sketched in Fig. ^. We briefly summarize the 
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Fig. 2. The phase diagram for an itinerant electronic system close to a 
quantum phase transition controled by a parameter b. The dashed lines 
refer to cross-overs. The region relevant for this paper is located to the left 
of be- 

results by Millis relevant for this paper. To the left of the quantum critical 
point the phase diagram is dominated by three regions. On lowering the 
temperature one first enters a regime which is dominated by classical order 
parameter fiuctuations. Below the doted line which indicates the Ginzburg 



temperature classical non-Gaussian order parameter fluctuations occur. De- 
creasing the temperature further the behavior on the model depends on 
the symmetry class of the fluctuations and whether the fermions are fully 
or partly gapped. Paramagnon like theories are often used in a parameter 
regime close to a 3D transition which is intercepting the low temperature 
regime of the quantum phase transition. It is often assumed that the or- 
der parameter fluctuations are still gaussian. However in a regime where 
the electronic density of states is already substantially suppressed like in the 
pseudogap regime this is a questionable assumption. The purpose of this pa- 
per is draw attention to the fact that the non-Gaussian fluctuations modify 
the electronic properties drastically compared to the Gaussian fluctuations. 



3. A QUASI-ONE-DIMENSIONAL PROBLEM 

To illustrate the effect of the non-gaussian fluctuations we discuss a 
model system for which most properties can be calculated without further 
approximation. Contrary to the assumption in previous work 00, it is not 
sufficient to describe the order parameter fluctuations by the variance and 
correlation length only, but that one needs to consider higher moments of the 
order parameter correlator. This becomes intuitively clear if one considers 
the situation depicted in Fig. |3|. If the order parameter varies smoothly. 
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Fig. 3. Typical order parameter configuration in the non-Gaussian regime. 
The electronic wave function is localized in the regions where the order- 
parameter is small. The overlap between electronic wave function localized 
in different regions is exponentially small. 

as in the Gaussian regime, regions where the order parameter is suppressed 



are smeared out over the correlation length. The electronic wavefunction 
is spread out over a length comparable with the correlation length. The 
kinetic energy is low and consequently many states can be found at low 
energy even when the correlation length is large. On the other hand, if the 
order parameter is established and only suppressed over a length scale much 
shorter than the correlation length of the potential, as in the non-Gaussian 
regime, the electronic wavefunction decays over a distance vp/A and has 
a large kinetic energy. For the same correlation length and variance of the 
order parameter the electronic wavefunction is much stronger suppressed for 
non-Gaussian fluctuations. The overlap of parts of the wavefunction located 
at regions where the order parameter is exponentially small. In dimension 
d the electron wave function will be constrained to a d — 1 dimensional 
manifold. If the dimension is greater than one the regions where the order 
parameter is small form a "network" in which the electrons can move and 
one has to consider "quantum-percolation". 

To be specific we consider here only the most simple case namely a 
one-dimensional system close to a CDW transition. The dispersion of the 
electrons close to the Fermi energy can be assumed to be linear. The Hamil- 
tonian of the electrons has the form: 



where the operators and create left and right moving electrons respec- 
tively. The classical order parameter field A{x) is determined by a Ginzburg 
Landau action given below, vp is the Fermi velocity. 

Next we consider the order parameter fluctuations. For commensurate 
fluctuations, the low energy electronic density of states is dominated by the 
Dyson singularity which only exists in one dimension I13'EJ. For the more 
general case, the order parameter fluctuations are complex and the Dyson 
singularity is absent. Therefore we will restrict our discussion to complex 
order parameters. After integrating out the electrons the classical complex 
order parameter fluctuations, A{x) are described by the Ginzburg-Landau 
action: 



Close to the mean field phase transition a varies linearly with the temper- 
ature a = a'{T/T^^ — 1), whereas P and Eire nearly temperature inde- 
pendent. In principle, the coefficients « and /? have to be determined 
self consistently from the electronic properties. The correlator for the or- 
der parameter fluctuations always decays exponentially, < A(x)A(0) >=< 
> exp(— since the ID system is disordered above the 3D ordering 



H = -ivF{R^d^R - L^d^L) + A{x)R^L + A*(a; 



)Lti?, 



(6) 





temperature, T^^ . Nevertheless the action Eq. (0) has two different regimes: 
if a{T) is positive and large, the order parameter fluctuations are centered 
around zero and basically Gaussian. For a(T) negative and large the am- 
plitude of the order parameter is given by \/< > and only the phase 
fluctuations play a role. 

The density of states of the electrons in an order parameter potential 
given by Eq. (j^ can be calculated exactly numericallv using the tranfer- 
matrix formalism for the order-parameter fluctuations^. The Schrodinger 
equation for the electrons right moving and left moving electrons, with the 
wavefunction u{x) and v{x) respectively, reads: 

— ivFdxu{x) + A {x)v{x) = eu{x) (8) 
+ivFdxv{x) + A*{x)u{x) = eu{x) (9) 

where e is the single-particle energy. The logarithms of u{x) and v(x) obey 
the following equation of motion: 

—ivFdxlog{u(x)) = e — A {x)u{x)/v{x) 
+ivFdxlog{v{x)) = e — A*{x)v{x)/u{x) 

Where one has to keep in mind choosing the correct branch for log(2;). 
Adding the two equation we obtain an equation for the ratio r{x) = u{x) /v{x): 

- ivF log(r(j;)) = 2 [e - A{x)r{x) - A*(j;)r~i(x)] (10) 

With the ansatz r(x) = exp(i(^(j;)) and A(x) = |A(a;)|e** we finally obain 
the equation of motion for the '^{x): 

dx^ix) = 2 [e - |A(x)| cos{^{x) + $(x))] (11) 

Note that ip{x) is the scattering phase of the wavefunction. Surprisingly the 
phase information alone is quite sufficient to obtain all the thermodynamic 
information. The integrated density of states, AA(e), can be related to the 
phase bya 

A^(^) = lim ^ (12) 

The density of states can then be obtained by numerically differentiating 
M{e) with respect to e. Calculating (p{x) for a particular configuration of 
the the order-parameter field amounts to integrating a one-dimensional dif- 
ferential equation with a random potential A{x). 

The remaining problem is to generate a typical configuration of A(x) 
and sum over all configurations with the correct statistical weight. In prin- 
ciple this can be done using a classical Monte-Carlo simulation. It turns out 
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Fig. 4. The electronic density of states as a function of temperature. The 
temperatures are 0.1, 0.2, 0.3 ... 1.0 x . The Ginzburg temperature 
is set to Tg = O.IT^^. The square boxes are the result of the Gaussian 
calculation in the limit correlation length going to infinity. 



that there is much more efficent way for this problem since eigenfunctions of 
the the transfer matrix for the free energy function Eq. |^ can be calculated 
easily. Using the transfer matrix formalism this problem can be mapped to 
a Langevin equation for A(x): 

^ = expH<I>)/(|A|)+r? (13) 
where / is the "guiding functions" of the random walk and rj is complex 



Gaussian white noise. If /(| A|) is a linear function then Eq. 13 just describes 
an Orenstein-Uhlenbeck process with a Gaussian spatial correlation. It has 
been shown that / can be obtained from the transfer-matrix of the Eq. |7[ 

The exact result for the density of states of the electrons as a function 
of temperature is shown in Fig. ^ and compared with the result of the cal- 
culation assuming Gaussian fluctuations with a very long correlation length 
0. Whereas the Gaussian model still allows many states in the mean- field 
gap the non-Gaussian fluctuations remove the states much more effectively. 
The artefact in the Gaussian model that the density of states behaves like 
A^(e) ~ even when the system is ordered is due to the fact that in the 



Gaussian model the most probable value for the order parameter is |A| = 
which is not true in the non-Gaussian model. In the Gaussian theories the 
surpression of the electronic states is due to large amplitude of the order 
parameter fluctuations which are strongly overestimated if one would like 
to obtain the same loss of spectral weight as in the non-Gaussian model. 
Since the paramagnon theories are Gaussian the naive application to the 
pseudogap regime is highly questionableElS. 



4. QUASIPARTICLES IN A TWO-DIMENSIONAL 
SUPERCONDUCTOR 

In this section we apply the methods developed in the previous chapter 
to the thermodynamic properties of a quasi-particle in a two-dimensional su- 
perconductor. It has been claimed that the thermal and transport properties 
of a high temperature superconductor above the transition are mostly de- 
termined by the phase-fluctuations of the preformed superconducting pairs. 
We will apply the quasi-classical approximation which is not really justified 
if the length scale of the phase fluctuations of the order parameter is com- 
parable to the lattice spacing. But in that case we will recover the "normal" 
quasi-particles anyways. 

It has been noted by SchopohlEl that the quasi-classical Eilenberger 
equation of a superconductor can be mapped to the linearized Bogoliubov- 
DeGennes (BdG) equation. If we identify u and v with particle and hole 
excitation in a superconductor Eq. ^ can be viewed as a BdG equation. 



With the phase equation, Eq. 11 it then is feasible to calculate the density 
of states of an inhomogeneus superconductor relatively easily. We have ap- 
plied this method to calculate the density of states of a quasiparticle put 
into a two-dimensional superconductor. We assume that the order parame- 
ter amplitude is already fixed. The phase fluctuations are then described by 
an xy-model. The classical xy-model can be simulated efficiently even close 
to the Bereshinskii-Kosterlitz-Thouless (BKT) transition using the Wolff 
algorithmii. We used a 128 x 128 grid for the xy - model. For each ther- 
mal configuration the phase equation is solved and the density of states 
is calculated. Finally we assume a d^^^yi gap - which is the relevant gap 
structure for the high temperature superconductors. The result of the cal- 
culation is shown in Fig. The temperature is varies from much below the 
BKT transition to very high temperatures. The suprising fact is that the 
BKT transition has very little effect on the equilibrium electronic degrees 
of freedom. At very high temperatures on does recover most of the density 
of states although the order parameter is already established. At very low 




Fig. 5. Density of states of a quasiparticle in a two-dimensional dj.-2_y2 
superconductor. The inverse temperature (3 = 1/T is varied from 0.1 to 2.0 
in steps of 0.1. The BKT transition is located approximately at J (5c ~ 0.7, 
where J is the coupling of the phases. The energy is measured in terms of 
the amplitude of the order parameter. 



sity of states for a (1^2 _y2 su- 



it would be interesting to 



temperatures one recovers the usual d-wave 
per conductor, which is linear at low energies 
apply the quasi-classical approach to the transport properties in the phase- 
fluctuation regime. However for transport vortices have a highly nontrivial 
effect on the electrons. 



5. CONCLUSIONS 

Paramagnon like theories have been used extensively to study electronic 
properties of itinerant electronic systems close to a phase transition. These 
theories amount to using a Gaussian approximation for the incipient order 
parameter field. In this paper we have discussed the effect of non-Gaussian 
order parameter fluctuations on the electronic properties. The non-Gaussian 
order parameter fluctuations have been shown to be relevant for the under- 
standing of a two-dimensional itinerant system close to a quantum phase 
transition. For a one-dimensional system close to a charge-density-wave 
transition the electronic properties can be calculated exactly. The suppres- 
sion of the density of states is much more dramatic than in Gaussian models. 



Finally we have investigated a quasiparticle in a two-dimensional supercon- 
ductor and mapped the problem to a one-dimensional BdG equation. 

Dedication: This paper is dedicated to Peter Wolfle on the occasion of 
his 60th birthday. 
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